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Abstract 

We derive the evolution equations of parton distribution functions appro- 
priate in different kinematic regions in a unified and simple way using the 
resummation technique. They include the Gribov-Lipatov-Altarelli-Parisi 
equation for large momentum transfer Q, the Balitskii-Fadin-Kuraev-Lipatov 
equation for a small Bjorken variable x, and the Ciafaloni-Catani-Fiorani- 
Marchesini equation which unifies the above two equations. We propose 
a modified Balitskii-Fadin-Kuraev-Lipatov equation from the view point of 
the resummation. This modified version possesses an intrinsic Q dependence 
and its predictions for the structure function Fi(x, Q 2 ) are consistent with 
the HERA data. 
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I. INTRODUCTION 



Perturbative QCD (PQCD), as a gauge field theory, involves large log- 
arithms from radiative corrections at each order of the coupling constant 
a s , such as InQ in the kinematic region with a large momentum transfer 
Q and ln(l/x) in the region with a small Bjorken variable x. These loga- 
rithms, spoiling the perturbative expansion, must be organized. To organize 
the logarithmic corrections to a parton distribution function, the various 
evolution equations have been derived. For example, the Gribov-Lipatov- 
Altarelli-Parisi (GLAP) equation |1J sums single logarithms InQ for a large 
x to all orders, the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation || sums 
ln(l/x) for a small x, and the Ciafaloni-Catani-Fiorani-Marchesini (CCFM) 
equation ||, appropriate for both large and small x, unifies the above two 
equations. 

The conventional derivation of the evolution equations usually requires 
complicated diagrammatic analyses. The idea is to locate the region of the 
loop momenta flowing through the rungs (radiative gluons) of a ladder dia- 
gram, in which leading logarithmic corrections are produced. For the GLAP, 
BFKL, and CCFM equations, the important regions are those with the strong 
transverse momentum ordering, the strong rapidity ordering, and the strong 
angular ordering, respectively. Summing the ladder diagrams with the differ- 
ent kinematic orderings to all orders, we obtain the evolution equations that 
organize the leading logarithms. The derivation will become very difficult, 
if one intends to improve the accuracy of the evolutions to next-to-leading 
logarithms. 

In this paper we shall propose an alternative approach to the all-order 
summation of the various large logarithms. This approach is based on the 
resummation technique, which was developed originally for the organization 
of double logarithms In 2 Q [Q. Recently, we applied this technique to some 
hard QCD processes, such as deep inelastic scattering, Drell-Yan production, 
and inclusive heavy meson decays, and demonstrated how to resum the dou- 
ble logarithms contained in parton distribution functions into Sudakov form 
factors || . It has been shown [^, [(J that the resummation technique is equiva- 
lent to the Wilson-loop formalism for the summation of soft logarithms. Here 
we shall further show that it can also deal with the single-logarithm cases, 
i.e. the evolution equations mentioned above. Therefore, the resummation 
technique indeed has a wide application in PQCD. 
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The procedures of the resummation are briefly summarized below. The 
derivative of a parton distribution function with respect to Q or x is first 
related to a new function, which contains a special gluon vertex. Express- 
ing this new function as the convolution of the subdiagram involving the 
special vertex with the original parton distribution function, we arrive at 
the evolution equation. The subdiagram, after factorized out of the new 
function, is identified as the corresponding kernel in the different kinematic 
limits. We show explicitly that the lowest-order subdiagram gives the kernel 
for the leading-logarithm summation exactly. Using our approach, it is not 
necessary to go into the detailed analysis of the orderings of radiative gluons 
in a ladder diagrams, since all the possible orderings have been included in 
the new function. Hence, the resummation technique provides a unified and 
simple derivation of the known evolution equations in the literature. 

In the region with both large Q and small x many gluons are radiated 
in scattering processes with small spatial separation among them, and a 
new effect from the annihilation of two gluons into one gluon becomes im- 
portant. Taking into account this effect, a nonlinear evolution equation, the 
Gribov-Levin-Ryskin (GLR) equation |7j is obtained. Using the resummation 
technique, the annihilation effect is introduced through next-to-leading-twist 
contributions to the subdiagram containing the special vertex, and the GLR 
equation can also be derived easily M. In the present work we shall not 
address this subject, because it is not very relevant. 

It is known that the BFKL equation is independent of Q, and thus its 
predictions are insensitive to the variation of Q. However, the recent HERA 
data |§ for the structure function F 2 (x, Q 2 ) of deep inelastic scattering ex- 
hibit a stronger Q dependence. It is plausible that the experiments have not 
yet explored the region with a low enough x such that the BFKL equation is 
applicable. Hence, to explain the data, the Q dependent GLAP equation is 
combined in some way. One can employ either the improved splitting func- 
tion, which embeds the BFKL summation, in the GLAP equation [PJ, or the 
CCFM equation directly |IT] , which unifies the BFKL and GLAP equations. 
In this paper we shall not rely on the GLAP summation of InQ, but propose 
a modified BFKL equation from the viewpoint of the resummation, which 
contains an intrinsic Q dependence. It will be shown that the gluon distribu- 
tion function derived from the modified BFKL equation gives the predictions 
of F 2 (x,Q 2 ), that are well consistent with the data. 

We derive the GLAP, BFKL and CCFM equations in Sections II, III, and 
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IV, respectively, by means of the resummation technique. We explain how 
the subdiagram containing the special vertex reduces to the corresponding 
evolution kernels in the different kinematic regions. The modified BFKL 
equation along with its analytical solution of the gluon distribution function 
are presented in Section V. We then evaluate the structure function F 2 (x, Q 2 ) 



using the high-energy p^-factorization theorem ||12|| , and compare them with 
the HERA data. Section VI is the conclusion. 



II. THE GLAP EQUATION 

We derive the GLAP equation in this section. Consider deep inelastic 
scattering (DIS) of a hadron with a light-like momentum p = p + 5^ + in the 
large x limit, where x = —q 2 /(2p ■ q) = Q 2 /(2p ■ q) is the Bjorken variable 
with q the momentum transfer to the hadron through a virtual photon. It 
is known that the collinear region with loop momenta of radiative gluons 
parallel to p is important, from which large logarithms arise ||. The other 
important regions are soft, with loop momenta much smaller than Q, and 
hard, with loop momenta being of order Q. In the covariant gauge dA = the 
scattered quark line the collinear gluons attach can be replaced by an eikonal 
line along an arbitrary direction n, n 2 / |. The associated Feynman rules 
are l/(n ■ I) for the eikonal propagator, and for an vertex on the eikonal 
line, I being the momentum flowing through it. With the eikonalization, the 
collinear gluons are factorized out of the hard scattering subamplitude and 
absorbed into a quark distribution funtion </>(£, p + , p) with £ the momentum 
fraction, and p a factorization or renormalization scale. The argument p + 
denotes the large logarithms ln(p + / p) appearing in 0, which will be organized 
below. 

The standard definition of <fi in the covariant gauge is then given by 

= / ^e-^-(p|g(y-)Vpe i /o ! '^-^)g(0)b) , (1) 

J Z7T Z 

as shown in Fig. 1(a), where 7 + is a Dirac matrix, and \p) denotes the in- 
coming hadron. The path-ordered exponential p e l \ dzn - A represents exactly 
the eikonal line collecting the collinear gluons. It is easy to confirm that 
this exponential generates the Feynman rules stated above. Though con- 
structed in this way contains an artificial dependence on n, the cross section 
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of the scattering process, being a physical quantity, should be n indepen- 
dent. It has been shown that the n dependence of is canceled by those of 
other involved subprocesses This vector n is, however, essential at the 
intermediate stage of the resummation. 

In the axial gauge n ■ A = the path-ordered exponential is equal to the 
identity, and is defined by Fig. 1(b), implying that the collinear gluons are 
decoupled from the eikonal line. The n dependence then goes into the gluon 
propagator, (-i/l 2 )N^ v (l), with 

N v» = gv» + n 2 —- - . 2 

n ■ I [n-iy 

The above decoupling can also be understood from the relation n^N^ = 
associated with the attachment of a radiative gluon to the eikonal line. We 
have shown that the resummation results derived in the axial and covariant 
gauges are the same |J, and thus we can work in either of them. In this 
section we adopt the axial gauge. 

The key step of the resummation is to obtain the derivative p + d<p/dp + . 
Because of the scale invariance of in the vector n as indicated by Eq. ([!]) 
or (0), must depend on p and n through the ratio (p ■ n) 2 /n 2 . Hence, we 
have the chain rule relating p + d/dp + to d/dn: 

, d , n 2 d 

P 77? = v a —(j) , (3) 

dp + v ■ n dn a 

with v a = 5 a+ a vector along p. In the axial gauge d/dn applies only to the 
gluon propagator, giving 

d -N^ = l —(FN au + FN^) . (4) 



dn a n ■ I 

The loop momentum Z M (l u ) flowing through the differentiated gluon line 
contracts with the vertex the gluon attaches, which is then replaced by a 
special vertex 



n 2 v r 



(5) 
v ■ nn ■ I 

This special vertex can be simply read off the combination of Eqs. (||]) and 
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The contraction of / M (l v ) hints the application of the Ward identity, from 
which it is found that the special vertex moves to the outer end of the valence 
quark line. For a given x, we obtain the formula, 

P + j^<l>{x,P + ,v) = 20(x,p + ,/i) , (6) 

shown in Fig. 2(a), where 4>, the new function mentioned in the Introduction, 
contains one special vertex represented by a square. The coefficient 2 comes 
from the equality of the new functions with the special vertex on either of 
the two valence quark lines. Note that Eq. (||) is an exact consequence of the 
Ward identity without specifying the ordering of radiative gluons. All the 
possible orderings are embedded in Eq. (|). 

We show that Eq. (|6|) leads to the GLAP equation, when p + , or equiva- 
lently Q in the center-of-mass frame of the virtual photon and the incoming 
quark, is large. The important regions of the loop momentum flowing through 
the special vertex are soft and hard, since the vector n does not lie on the 
light cone, and the collinear enhancements are suppressed. In the important 
soft and hard regions <p can be factorized into the convolution of the subdi- 
agram containing the special vertex with the original distribution function 



4>{x, p + , fx) = f d£[K(x, f , p + , fi) + G(x, f , p + , p + , //) . (7) 

The function K, absorbing the soft divergences of the sub diagram, corre- 
sponds to Fig. 2(b), where the eikonal approximation for the valence quark 
propagator has been made. The function G, absorbing the ultraviolet diver- 
gences, corresponds to Fig. 2(c), where the subtraction of the second diagram 
ensures that the involved loop momentum is of order Q. 

We identify the factor K+G as the GLAP kernel. According to Figs. 2(b) 
and 2(c), K and G are written as 



K = ig 2 C F Li e J 



dr H V..V- 



(2vr) 4 - e v ■ I 



+ 2m5(l 2 )5(t - x - l + /p + ) 



-8K , (8) 

r d 4 ~ e l 
G = -ig 2 C F ^J 



(2tt 



,4-e ^ 



{t,p-l) 2lv v-l 



__ S (C -x)-SG, (9) 
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where Cf = 4/3 is the color factor. The 8 functions come from the final 
state cut, and SK and SG are additive counterterms. The virtual and real 
gluon emissions correspond to the first and second terms in the integral of K, 
respectively. The second term in the integral of G is the soft subtraction. We 
emphasize that the factorization formula in Eq. (0) is not an exact relation, 
but holds only up to the leading logarithms lnp + . As explained later, the 
approximation equivalent to the strong transverse momentum ordering in the 
conventional approach has been applied to Eq. fl8|). 
A straightforward calculation gives 



& = ~Z 



1 + i n ^ 5 (i_ x /£) 



G = -^M Cf I^S(1- X /Q, (10 ) 

< /i 

where constants of order unity have been dropped, and v = J (v ■ n) 2 /\n 2 \ is 
the gauge factor. Equation (|T0|) confirms our argument that depends on p 
and n through the ratio (p ■ n) 2 /n 2 = (up + ) 2 . In the considered region with 
x — > 1 the logarithm ln(^z/p + //i) in G can be replaced by ln(i/p + //i). 
We then treat K and G by renormalization group (RG) methods: 

fi^-K = -X K = -n-fj-G . (11) 
Cl/i a/i 

The anomalous dimension of K is defined by Xk = —fidSK/dn, whose explicit 
expression is not essential here. When solving Eq. (|TT|), we allow the variable 
fi to evolve from the scale of K to the scale of G. The RG solution of K + G 
is given by 



K(x,^p + ,fi) + G(x,^,p + ,fx) = K(x,£,p + ,p + ) + G(x,£,p + ,p + ) 

fP+ ^\ ( (-W 

- t Cf fi;' (12) 

It is obvious that the source of double logarithms, i. e. the integral cantaining 
Xk, vanishes. 
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A remark is in order. The function 5(£ — x — l + /p + ) in Eq. (|j) will 
be replaced by 5(£ — x)exp(ilr • b), b being the conjugate variable of the 
transverse momentum carried by the valence quark, if the transverse degrees 
of freedom are considered ||. The integration over £ can thus be performed 
trivially, and the convolution form in Eq. (|7|) is simplified to a multiplication 
form, 

<j>(x, 6, p + , fi) = [K(x, b, fi) + G(x, p + , n)]<p(x, b, p + , fi) . (13) 

This simplication will be implemented in Sec. V to solve the modified BFKL 
equation. In this case the scale 1/6, instead of p + , serves as the infrared 
cutoff of K. Consequently, 1/b is substituted for the lower bound of p, in 
Eq. fljjP. Double logarithms then exist, implying that the soft logarithms in 
4> do not cancel completely. Therefore, the resummation technique can deal 
with both the single-logarithm and double-logarithm problems. 
Inserting Eq. QT2j) into (|7|) and solving (^), we obtain 

<t>(x, fj.) = <f>(x, A, n) + [" d ±^E±C F f f —J—M, ft, , (14) 

JA fi 71 Jx £ (1 — X/Q + 

with A an arbitrary cutoff. We have defined (p(x,fi) = <p(x, [/,, ji), which does 
not contain large logarithms. According to the formalism in |J, the upper 
bound of fl should be set to the hadron momentum. However, it is equivalent 
to set it to /i here. It is known that a physical quantity such as the DIS cross 
section is \x independent. When applying the RG analysis to the factorization 
formula for the cross section, which is a convolution of the hard scattering 
subamplitude with 0, \i evolves to the characteristic scale of the hard part 
at last. This hard scale is exactly the large hadron momentum. 

Differentiating Eq. ([TJ|) with respect to fi, and substituting the RG equa- 
tion nd(f)(x(£), A(/i), n)/dfM = — 2A g 0(a;(£), A(/i), /i), X q = —a s /7i being the 
quark anomalous dimension in the axial gauge, we have 

d ,. . Otsifl) n f 1 d£ 2 , / x,/ x ft i-\ 

/i— (p(x, fi) = C F / —- t— 0(f , fi) - \ q {fi)<p{x, fi) . (15) 

dfi 7i Jx 4 (1 - x/E > ) + 

The above equation can be reexpressed as 

Q-^, Q) = ^ t ^P(x/mt, Q) , (16) 

dU 7T Jx t 
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where p, has been set to Q and the kernel P is 



P(x) = C f 



x). 



(17) 



It is trivial to identify P as the splitting function P qq in the limit x — > 1, 



Pqq(x) 



C f 



1 



1+X2 L 3 AM ^ 



(18) 



Hence, Eq. (y) leads to the GLAP equation (jig). 

In fact, only the terms of P 9g , which are singular at x — *■ 1, can be repro- 
duced in the resummation formalism. We emphasize that the factorization of 
the subdiagram containing the special vertex makes sense only in the leading 
soft and hard regions, and thus the functions K and G collect only the most 
important contributions to the splitting function. This conclusion applies to 
other cases that involve a soft approximation, such as the CCFM equation in 
the conventional derivation [[|, where the finite part of the relevant splitting 
function was also missing and put in by hand at last. 



III. THE BFKL EQUATION 

In this section we demonstrate that the resummation technique reduces 
to the BFKL equation for the gluon distribution function in the small x 
region. It is convenient to adopt the covariant gauge OA = 0, under which 
the unintegrated gluon distribution function F(x,px) is defined by Fig. 1(a) 
with the valence partons being gluons. The path-ordered exponential in 
Eq. (HD, and thus the eikonal line in the direction n, which collects radiative 
gluons, appear. F(x,pr) describes the probability of a gluon carrying a 
longitudinal momentum fraction x and a transverse momentum p?. We have 
made explicit that F depends on p?, instead of the large hadron momentum 
p + appearing in the GLAP case. It will be shown later that p + disappears as 
x — > 0, and thus the pt dependence is not negligible. For a unified treatment 
of the large-x and small-x summation, the p + dependence should be included, 
and the resummation technique reduces to the CCFM equation discussed in 
Sec. IV. 

To sum large logarithms ln(l/x), the strong rapidity ordering of radiative 
gluons in a ladder diagram is assumed in the conventional approach. In 
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the resummation formalism the BFKL equation can be derived simply by 
reinterpreting the derivative with respect to p + in Eq. (^) and by modifying 
the expression of K in Eq. (g). Though F does not depend on p + explicitly, 
it can vary with p + through the momentum fraction implicitly, which is 
proportional to (p + )~ 1 for a fixed parton momentum. For a similar reason, 
F depends on the ratio (p ■ n) 2 /n 2 , and thus Eq. @ holds. In the covariant 
gauge the operator d/dn applies to the Feynman rules for the eikonal line, 
giving 

*J0L * U _ v*) . (19) 

dn a n ■ I n ■ I \ n ■ I J 

Combining Eqs. and (IS), we find that the differentiation with respect to 



p + generates a special vertex on the eikonal line 

n 2 ( v ■ I 



n. 



v • n \n ■ I 



n a -v a ) . (20) 



The derivative of F is then expressed as 

p + -f—F(x,p T ) = -x-^-F(x,p T ) = AF(x,pt) , (21) 
dp + dx 

described by Fig. 3(a), where the new function F contains one special vertex 
denoted by the symbol x. Note that the coefficient 4 in front of F is twice 
of the corresponding coefficient in the GLAP case. Since the gluon interacts 
with the virtual photon through a quark box, two quark lines, and thus 
two eikonal lines after factorizing out the gluon distribution function, are 
adjacent to the parton vertex. Hence, there is one more attachment of the 
special vertex to the eikonal line on each side of the final state cut. It is 
trivial to show that these attachments give the same results. Starting with 
Eq. (0), we need not to go into the detailed analysis of the rapidity ordering 
of radiative gluons, since all the possible orderings have resided in it. 

The leading regions of the loop momentum I flowing through the spe- 
cial vertex are also soft and hard. If I is collinear, the first term v ■ I in 
n a vanishes, and the second term contracted with a vertex in the 

distribution function which is dominated by momenta parallel to p, gives a 
small contribution. The soft divergences of the subdiagram are collected by 
Fig. 3(b), and the ultraviolet divergences by Fig. 3(c). We employ the re- 
lation fabctbtc — (i/2)Nt a for the color structure, t being the color matrices 
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and N = 3 being the number of colors. Absorbing t a into the parton vertex, 
Fig. 3(b) leads to 

dH Y^ vX h v 



soft 



>,Pt) 



2 y J (2tt) 4 (2v ■ l)n ■ I 
0(P 2 T - l 2 T ) 



27t^(/ 2 )F(x,| P t + 1t|) 



P 



-F(x,p T ) 



where the triple-gluon vertex for vanishing I is given by 



(22) 



(23) 



The first term in the above integral corresponds to the real gluon emission, 
where F(x, \p? + It I) indicates that the parton coming out of the hadron 
carries a transverse momentum + 1^ in order to radiate a real gluon of 
momentum 1^. The second term corresponds to the virtual gluon emission, 
where the 9 function sets the upper bound of It to pt to ensure a soft 
momentum flow. Therefore, it is not necessary to introduce a renormalization 
scale /i into F. There is not the £ integration in Eq. (^), since we have 
included the transverse degrees of freedom of the parton as a soft regulator, 
and the final state cut associated with the real gluon can be approximated 
by <5(£ — x — l + /p + ) ~ <5(£ — x) as stated in the previous section. 

It can be easily shown that v x in Eq. (|23|) , contracted with a vertex in 
the quark box diagram, leads to a contribution smaller by a power 1/s with 
s = (p + q) 2 , compared to the contribution from the last term v v . Similarly, 
the term t> M is also contracted with a vertex in the box diagram through 
the metric tensor associated with the gluon distribution function. Hence, 
we drop the first two terms of T fluX , and absorb g Xfl into F . Evaluating the 
integral straightforwardly, Eq. (^2]) reduces to 



o;.s 



soft 



'x,p T ) = -j 



irl T 



F(x,\p T + l T \)-9(p 2 T 



l T )F(x,p T ) 



(24) 



with a s = Na s /7i. 

It was argued that when the fractional momentum of a parton vanishes, 
the associated collinear enhancements are suppressed ||. The vanishing of 
the contribution from the first diagram of Fig. 3(c) 



T 



d 2 i 



T 



7l 



In-* 



l\ + (xp + v) 
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xp+v , Jl T + (xp + u) 2 - xp + u 
In 



2 [/2 + (xp+J/) 2 ] 3 / 2 ^2 + ( xp+I/ )2 + xp + v 



(25) 



at x — ► 0, reflects this argument. Hence, F does not acquire a dependence on 
the large scale p + , and the transverse degrees of freedom must be taken into 
account, differing from the GLAP case for a large x. This is the basic idea 
of the so-called high-energy p^-factorization theorem [12]. Therefore, the 
introduction of px and the disappearence of p + are built in the resummation 
technique naturally. It is obvious that our formalism is applicable to the 
distribution functions constructed according to the collinear factorization 
(the GLAP case) and according to the p^-factorization (the BFKL case). 

Neglecting along with its soft subtraction (the second diagram in 
Fig. 3(c)), that is, adopting F = F so ft, Eq. ( PT| ) becomes 



F{x,p T ) _ f d% 

- a a 



J ^ [F(x, \p T + l T \) - 6(p T - l 2 T )F(x,p T )} , (26) 



d\\i{l/x) J ■kIj, 

which is exactly the BFKL equation. It is then understood that the subdia- 
gram containing the special vertex plays the role of the BFKL kernel. 

In summary, the BFKL equation is appropriate for the multi-Regge re- 
gion, where the transverse momenta flowing through the rungs of a ladder 
diagram are of the same order, i.e. It ~ Pt- Hence, the loop momentum 
It flowing through the parton distribution function is not negligible, and the 
final state cut <5(£ — x — 1 + /p + ) can be approximated by 5(£— x) as in Eq. (|22|). 
While the GLAP equation is appropriate for the transverse momemtum or- 
dered region, in which we have It Pt, i-e- F(x, |pr + 1t|) ~ F(x,pt) for 
the real gluon emission. The Pt dependence of the parton distribution func- 
tion then decouples, and can be integrated out from both sides of Eq. ( |26| ) . 
Hence, a parton distribution function in the GLAP equation does not involve 
the transverse degrees of freedom. In this case the loop momentum Z + should 
be maintained in the final state cut as in Eq. (§), which then leads to the 
splitting function. If not, the right-hand side of Eq. (Q) will be identical to 
zero. Therefore, the subdiagram containing the special vertex reduces to the 
corresponding evolution kernels in the different kinematic regions. 



IV. THE CCFM EQUATION 
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Based on the discussion in the previous two sections, it is not difficult to 
demonstrate that Eq. reduces to the CCFM equation ||, which is ap- 
propriate for both large x and small x. It embodies the GLAP equation and 
the BFKL equation, and depends on the longitudinal momentum p + and the 
transverse momentum px of a parton at the same time. For the conventional 
derivation of the CCFM equation, assuming the angular ordering of radia- 
tive gluons in a ladder diagram, refer to ||. By means of the resummation 
technique, the complicated diagrammatic analysis can be avoided, and the 
physical meaning of each factor in the CCFM equation is clearer. 

Consider Eq. fl2~l~|) but with the unintegrated gluon distribution function 
F depending on px and p + , 

P+ dp+ F< ^ X} PT} P+ ^ = 4 ^ X ' PT ' P+ ^ ' ^ 
which manifests the attempt to unify the GLAP and BFKL equations. It 
is not necessary to introduce a renormalization scale /i here, because the 
transverse degrees of freedom will not be integrated out. Again, the new 
function F involves one special vertex on the eikonal line in the direction n. 

If following the standard procedures of the resummation, we should fac- 
torize the subdiagram containing the special vertex by absorbing its soft and 
hard contributions into the functions K and G, respectively. The p + depen- 
dence of F then comes from G, which collects the virtual gluon emissions. 
This idea leads to a new unified evolution equation, which will be studied 
elsewhere. To reproduce the CCFM equation, however, the inclusion of the 
virtual corrections must be performed in a different way: We organize the 
virtual gluons embedded in K, instead of those in G. Hence, the subdiagram 
is factorized into Fig. 4(a), where the two jet functions J group all the pos- 
sible virtual corrections, and the real gluon between them is soft. It will be 
shown below that this subdiagram gives the CCFM kernel. 

First, we resum the double logarithms contained in J by considering its 
derivative 

P + ^ j (Pt,P + ) = 2J(p T ,p + ) 

= 2[Kj(p T ^) + Gj(p + lf i)}J(p T ,p + ) . (28) 

At lowest order the function Kj comes from the first diagram of Fig. 3(b), 
and Gj from Fig. 3(c). The coefficient 2 counts the two eikonal lines adja- 
cent to the parton vertex. The relation between Kj + Gj and J is simply 
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multiplicative, since J groups only virtual gluons. We have set the infrared 
cutoff of Kj to pt, as indicated by its argument. This cutoff is necessary 
here due to the lack of the corresponding real gluon emission, which serves as 
a soft regulator. The one-loop Kj can be easily obtained by working out the 
second integral in Eq. ( p4[) without the 9 function. The anomalous dimension 
of Kj is then found to be jj = a s /2. The function Gj can also be computed, 
but its explicit expression is not important. The standard RG analysis gives 



fP + dfi 

Kj(p T , fj) + Gj(p + , n) = - / — jj(a s (p)) , 



with the initial conditions Kj(p T ,p T ) = Gj(p + ,p + ) = 
neglected the constants of order unity in Kj and Gj. 
Substituting Eq. ( |29"D into (|28|), we solve for 

J(p T ,Q)=A(Q,p T )jW , 



with the double-logarithm exponential 



A{Q,p T ) = exp 



-a. 



dp + fP + dfi 



PT P JPT 



(29) 

0. Of course, we have 
(30) 

(31) 



We have chosen the upper bound of p + as Q, and ignored the running of a s . 
The initial condition regarded as the tree- level gluon propagator, will 
appear in the integrand for the real gluon emission in Fig. 4(b). We split the 
above exponential into 



A(Q,p T ) = A 1 s /2 (Q,zq)A^ s (z,q,p T ) 



1/2. 



(32) 



with z = x/£ and q = iy/(l — z), where £ is the momentum fraction entering 
J from the bottom, and It is the transverse loop momentum carried by the 
real gluon. The so-called "Sudakov" exponential A s and the "non-Sudakov" 
exponentials A are given by 



&s(Q,zq) 



exp 
exp 



-2a K 



-at. 



Q dp + rp + dfl 

zq p^ JpT 
1 J " 2 '1-Pt/p dz' 



dp 1 



{zq)* p* JO 



Z' 
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exp 
exp 



-2a, 



-a s 



"i dp + fP + d[i 

PT P + JPT P 

pt/q dz' fP 2 T dp 



z' J(z'q) 2 p* 



(33) 



where the variable changes p, = (1 — z')p and p + = p for As, and p, = p 
and p + = z'q for A^s have been employed. The inserted scale zq reflects 
the combination of the rapidity ordering for the BFKL equation and the 
transverse momentum ordering for the GLAP equation || |TT|j . 

Picking up the last term of the triple-gluon vertex in Eq. (|23|), F is writ- 
ten, in terms of Fig. 4(b), as 



F(x,p T ,p^ 



o 

x5(e 



dt 



d 4 l v^h, 



j2mS(l 2 )A 2 (Q,p T ) 



x 



{2-kYv- In ■ I 
l+/p+)9(Q - zq)F(£, |p T + 1t|,P + ), (34) 



with / the gluon momentum. The propagator 1/v ■ I comes from the eikonal- 
ized tree-level on the right-hand side of Fig. 4(b). The left-hand side 
has been absorbed into F. The argument of F in the integrand is pt + It, 
because it is for the real gluon contribution. Basically, the above formula 
is similar to the real gluon part of the BFKL equation ( |22|) except for the 
exponential A 2 from the two jet functions J, and 6(£ — x — 1 + /p + ) associated 
with the final state cut, which is restored due to the inclusion of the p + de- 
pendence as in Eq. (|J). Through this 5 function, the GLAP splitting function 
is generated. The 9 function guarantees that the Sudakov exponential A s is 
meaningful. Hence, the expression of Eq. ( |34]) manifests the combination of 
the GLAP and BFKL features. 

Performing the integration over l + and l~, we obtain 



F(x,p T ,p A 



a t 
1 



d 2 h 



2n 2 {i-x)p +2 



7T [n + l 2 T + 2n-(i 
x9(Q-zq)F(Z,\p T + l T \,p + ) , 



x) 2 p +2 ] 2 



A 2 (g,p T ) 



(35) 



where we have assumed n = (n + ,n~,0) for convenience. Eq. ( |35D is then 
substituted into (|27|) to find the solution of F. We adopt the variable changes 
£ = x/z and It — (1 — z)q, and integrate Eq. Q3"5D from p + = to Q. 
To work out the p + integration, F(x/z, \p? + 1t|,P + ) is approximated by 
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F(x/z, |pt + It\,It)- This approximation is fine, if one does not intend to 
obtain the ^-function terms of the splitting function P gg , written as 



P99 ~ a $ 



1 



1 



+ - + z (l- z ) + ^-N--n f )5(l-z) 



1 — z z 

rif being the number of quark flavors. Eq. ( p5|) then becomes 

"l r d 2 q 



(36) 



F{x,pr,Q) = F^ + a s [ dz f -^-6{Q-zq)As(Q,zq)A NS {z,q,p T 

Jx J Tcq 



X- 



-F{x/z,\p T + {l-z)q\,l T ) , (37) 



z{l-z, 

where the nonperturbative initial condition F^ corresponds to the lower 
bound of p + , and the term suppressed by l/Q 2 in the integral has been 
dropped. 

Eq. ( j3~7|) can be rewritten as 



(38) 



F(x,pr,Q) = F®+[ dz [^6(Q-zq)A s (Q,zq)P(z,q,p T ) 

Jx J nq z 



xF(x/z, |p T + (1 - z)q\,l T ) 



with the splitting function 
P = a 



— h A NS (z,q,p T )- + z(l - z) 

1 — z z 



(39) 



To arrive at the above expression we have employed the identity l/(z(l — 
z)) = 1/(1 — z) + 1/z, and neglected the non-Sudakov form factor Ajvs in 
front of 1/(1 — z), because A^s has no pole as z — > 1 as shown in Eq. (|33|). 
Obviously Eq. (fEJ) is the CCFM equation §. The last term z(l - z) of P 
is put in by hand, which is hinted by Eq. (|36|). This term, finite at z — > 
and at z — > 1, can not be obtained in the conventional approach either |§ as 
stated at the end of Sec. II. 



V. A MODIFIED BFKL EQUATION 

We have mentioned in the Introduction that the recent HERA data of the 
DIS structure function F 2 (x,Q 2 ) can be explained by the CCFM equation 
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11|], in which the Sudakov exponential A$ collects the summation of InQ 
for a large x. In this section we propose a modified BFKL equation based on 
the resummation technique. This modified equation, much simpler than the 
CCFM equation, contains an intrinsic Q dependence, instead of that from 
the InQ summation. It will be shown that the resultant predictions for F 2 
match the HERA data. With this success, we demonstrate the power of the 
resummation technique. 

An alert reader may have noticed that to derive the BFKL equation, 
we must extend the loop momentum l + to infinity in Eq. ( pffi) . Strickly 
speaking, the real gluon emission in fact involves the distribution function 
F(x + l + /p + , |pt + 1t|) as part of the integrand, which is then approximated 
by F(x, |pt + 1t|) in the soft I region. Therefore, the behavior of F, vanishing 
at the momentum fraction equal to unity, should introduce an upper bound 
of Z + . To obtain a more reasonable BFKL kernel, we truncate l + at some 
scale, and a plausible choice of this scale is of order Q. We then derive a 
modified BFKL equation for a Q dependent gluon distribution function, 

dF(x,p T ,Q) _ [ d 2 l 
a 



J -p[F(x, |p T + 1t|, Q) - 0{Ql - l 2 T )F(x,p T , Q)\ 



d\n(l/x) J 

in the x — > limit, where the last term comes from the upper bound of / + . 
The gauge vector n has been chosen to render the coefficient of Q 2 equal to 
unity. It can be shown that our predictions are insensitive to this coefficient, 
as long as it is of order unity. Obviously, Eq. (fJ(|) approaches Eq. (^) in the 
Q —>■ oo limit. 

Another modification is that the loop momentum Ij- in the virtual gluon 
emission is truncated at an arbitrary scale Qq of order 1 GeV, instead of 
Pt as in the conventional BFKL equation. This modification is reasonable, 
since the virtual gluon contribution plays the role of a soft regulator for 
the real gluon emission only, and setting the cutoff to Qq serves the same 
purpose. Furthermore, the replacement of pr by Qq allows us to solve Eq. ([40|) 
analytically, and the solution of F maintains the essential BFKL features. It 
will be shown that the first term of the integral is responsible for the rise of 
F, and the last term acts to slower the rise. We then expect that F ascends 
faster at a larger Q, for which the effect of the last term is weaker. 
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We observe that the CCFM equation will be identical to the modified 
BFKL equation, if the exponential A 2 is "unfolded" into the lowest-order 
virtual gluon contribution, Q(Q—zq) is dropped, |p T +l T |,p + ) is replaced 
by F(x, \pt + 1t|,Q), and the £ integration is performed to produce the Q 
dependent term. Therefore, the Q dependence of the latter is not attributed 
to the all-order InQ summation in A, but to the boundary of the phase space 
for radiative corrections. In fact, this Q dependent term also appears in the 
CCFM equation, but was dropped in the derivation of Eq. (|37|). 

The Fourier transform of Eq. (|40|) leads to 



dF(x,b,Q) 
d\n(l/x) 



d 2 ln 



71 



a s (l/b) 
-S(b,Q)F{x,b,Q) , 



I'rp 



F(x,b,Q) 



(41) 



with 



S(b,Q) = 2a s (l/b)[\n(Q b)+ 1 -\n2 + K (Qb)} . (42) 

The Bessel function K comes from the last term of the integral, and 7 is the 
Euler constant. The argument of a s has been chosen as 1/6, since we work 
in the conjugate b space. 

Eq. fljip can be trivially solved to give 



F(x, b, Q) = F(x , b) exp[— S(b, Q) \n(x /x)\ 



(43) 



Xq being the initial momentum fraction below which F begins to evolve ac- 
cording to the BFKL summation. Transforming Eq. ( |4*3"D back to the momen- 
tum space, we derive the analytical solution to the modified BFKL equation, 



F(x,p T ,Q) 

and the gluon density xg by integrating Eq. 
d 2 p T 



bdbJ {p T b)F{x,b,Q) 



(44) 



over p T , 



xg(x,Q 



71 



F(x, p T} Q) = 2Q / dbJ x (Qb)F(x, 6, Q) . (45) 



In the above expressions Jo and J\ are the zeroth and first order Bessel 
functions, respectively. 
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To proceed with the numerical analysis, we assume a "flat" gluon distri- 
bution function [|ll], [13| 

F^(x,b) = 3N g (l-xfexp(-Q 2 b 2 /4) , (46) 

for x > xq, N g being a normalization constant, which is the Fourier transform 
of 

F { °\x,Pt) = ^N g (l - xfeM-Pr/Ql) (47) 

in momentum space. The initial condition F(xq, b) in Eq.(f43|) is then equated 
to F(x ,b) = F(°)(x ,&). We set Q = 1 GeV and x = 0.1 arbitrarily. It 
can be shown that the predictions vary only slightly for other choices of the 
parameters Qq and Xq of the same order. N g will be determined by the 
data of the structure function F 2 (x,Q 2 ) at a specific value of Q 2 , and then 
employed to make predictions for other values of Q 2 . 

An advantage of the b space is that the infrared sensitivity of the BFKL 
solution from the pt diffusion is moderated. The divergence of a s (l/b) from a 
large b is suppressed by the exponential e - *^ 2 / 4 in the initial condition P(°\ 
In the momentum space, however, the divergence of ck s (pt) from small px 
is not suppressed by exp(— p^/Ql) in F^°\ and thus the solution is sensitive 
to the infrared cutoff of pt- We have confirmed that our predictions of F 2 
almost remain the same for the cutoff of b at 2-4 GeV -1 . 

We compute the structure function F 2 , whose expression, according to 
the pT-factorization theorem, is given by 

F 2 (x,Q 2 ) = [ Pc fPZH( X /t,p T ,Q)F(t,p T ,Q) , (48) 

Jx t, JO 71 

p c being the upper bound of px which will be specified later. The hard 
scattering subamplitude H denotes the contribution from the quark box dia- 
grams, where both the incoming photon and gluon are off shell with q 2 = —Q 2 
and p 2 = —px, respectively. For simplicity, we consider only the contraction 
—g^Wnv in the calculation of H, and neglect the contribution iiomp^p u W llu , 
which is less important, being the DIS hadronic tensor. We also assume 
that a charm quark is massless, and that a bb quark pair is not involved in 
the box diagram, i.e. the active flavor number rif in the running coupling 
constant a s is equal to 4. Following these assumptions, we concentrate on 
the range of Q 2 between 8 and 20 GeV 2 . 
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A simple calculation gives 



H{Z,P T ,Q) = e2 q^ Z< 



, 2 + (l-z) 2 -2,(l-2z)^ + 2, 2 |: 



1 1 + J 1 - Az 2 p 2 T /Q 2 ) 
x - In " — — 2} (49) 

'l - Az 2 p 2 T /Q 2 l-Jl- 4z 2 p 2 T /Q 2 J 



with e q the electric charge of the quark q. Note that the terms in the braces 
approach the splitting function 

P qg (z) = \[z 2 + (1 - z) 2 ] (50) 

in the pt — > limit. To require a meaningful H, we modifiy the upper bound 
of pt in Eq. ( |48|) from p c = Q to 

Pc = min (q, J-Q^J . (51) 

We evaluate the integral in Eq. fl48|) straightforwardly for Q 2 = 15 GeV 2 , 
and then determine the normalization constant N g = 3.656 from the data 
fitting. When Q 2 varies, we adjust N g such that xg has a fixed normalization 
Jq 1 xgdx. It is found that N g changes only by about 5% in the considered 
range of Q 2 . F 2 for Q 2 = 8.5, 12, and 20 GeV 2 are then computed, and 
results along with the HERA data |§ are displayed in Fig. 5. It is obvious 
that our predictions agree with the data well. The curve has a steeper rise 
at a larger Q, which is the consequence of the Q dependent modified BFKL 
equation. The curves descend rapidly at x close to 0.1, since the contributions 
from other kinds of partons, such as the valence quarks, which are more 
important in this intermediate x region, are not included. For comparision, 
we also present the results from the conventional BFKL equation, which can 
be obtained simply by substituting l T + M 2 for the denominator l T + Q 2 
in Eq. (|40|), or equivalently, Mb for the argument Qb of the Bessel function 
K in Eq. ( fT2"| ) with an extremely large M = 10 3 -10 4 GeV. In this case the 
normalization constant determined from the best fit to the data for Q 2 = 15 
GeV 2 is N g = 2.908. It is found that the shape of the curves is almost 
independent of Q, and thus the match with the data is not very satisfactory. 
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At last, we show in Fig. 7 the behavior of the gluon density xg computed 
from Eq. ([45]). The rise of xg at small x is due to the term InQob hi the 
exponent S as mentioned before, which leads to an integrand proportional 
to 

(Q b)- 2aaHxo/x) . (52) 

The small b region then gives a huge contribution. When x approaches zero 
such that 2a s \n(xo/x) > 1, the integration over b diverges. However, the 
last term of S, Ko(Qb) oc — ln(Qb) in the b —> limit, which comes from the 
upper bound of the loop momentum l + , cancels the divergence. This is the 
reason l + must be truncated. On the other hand, the variation of xg with x 
for different Q can also be understood from the combination of ln(Q b) and 
K (Qb) in the b — > limit, written as, 

/ rf, \ -2a s ln(<2/Qo) 

(I) ' (53) 

The above expression indicates that the exponent A, characterizing the rise 
of xg ~ x~ x at small x, increases with Q. The values A ~ 0.36 for Q 2 = 8.5 
GeV 2 and A ~ 0.51 for Q 2 = 20 GeV 2 are deduced from Fig. 6, which 
are consistent with that obtained in ]14| from a phenomenological fit to the 
HERA data (A « 0.3 for Q 2 = 4 GeV 2 ) and with that in |0| from solving 
the conventional BFKL equation numerically (A ~ 0.5 for a wide range of 

Q 2 ). 



VI. CONCLUSION 

In this paper we have shown that the resummation technique provides 
a unified and simple viewpoint to the organization of the various large log- 
arithms, and reduces to the GLAP equation, the BFKL equation, and the 
CCFM equation in different kinematic regions. The main idea is to relate 
the derivative of a parton distribution function to a new function involving a 
special vertex. The summation of the large logarithms is then embedded in 
the new function without resort to the complicated diagrammatic analyses. 
When expressing the new function as a factorization formula, we obtain the 
evolution equation, and the subdiagram containing the special vertex is ex- 
actly the corresponding kernel. By means of the resummation technique, the 
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derivation of the evolution equations is simpler. Furthermore, to improve 
the accuracy of the kernel to next-to-leading logarithms, we only need to 
evaluate the 0(a 2 s ) subdiagram. Such an evaluation can be performed in a 
straightforward way. The BFKL equation including the summation of the 
next-to- leading ln(l/x) will be published elsewhere. 

We have also calculated the DIS structure function F 2 (x,Q 2 ) using the 
analytical solution of the unintegrated gluon distribution function from the 
modified BFKL equation. Our predictions exhibit a stronger Q dependence, 
and are in a good agreement with the HERA data, compared to those from 
the conventional BFKL equation. Note that the Q dependence in the modi- 
fied BFKL equation is intrinsic, which arises from the boundary of the phase 
space for radiative corrections, instead of from the InQ summation in the 
CCFM equation. In this sense we argue that the current experiemnts may 
have explored the multi-Regge region with ln(l/x) ^> InQ. Certainly, this 
issue still needs to be clarified by further theoretical and experimental studies 
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Figure Captions 



FIG. 1. Definition of a parton distribution function in (a) the covariant 
gauge and in (b) the axial gauge. 

FIG. 2. (a) The derivative p + d<f)/dp + in the axial gauge, (b) The 0(a s ) 
function K. (c) The 0(a s ) function G. 

FIG. 3. (a) The derivative —xdF/dx in the covariant gauge, (b) The soft 
structure and (c) the ultraviolet structure of the 0(a s ) subdiagram contain- 
ing the special vertex. 

FIG. 4. (a) The subdiagram containing the special vertex for the CCFM 
equation, (b) The subdiagram for the CCFM equation after resumming the 
double logarithms in J. 

FIG. 5. The dependence of F2 on x derived from the modified BFKL equa- 
tion (solid lines) and from the conventional BFKL equation (dashed lines). 
The HERA data || are also shown. 

FIG. 6. The dependence of xg on x derived from the modified BFKL equa- 
tion for, from bottom to top, Q 2 = 8.5, 12, 15, and 20 GeV 2 . 
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